Role of the van Hove singularity in the Quantum Criticality of the Hubbard Model 
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A quantum critical point is found in the phase diagram of the two-dimensional Hubbard model 
[Vidhyadhiraja et al, Phys. Rev. Lett. 102, 206407 (2009)]. It is due to the vanishing of the critical 
temperature associated with a phase separation transition, and it separates the non-Fermi liquid 
region from the Fermi liquid. Near the quantum critical point, the pairing is enhanced since the 
real part of the bare d-wave pairing susceptibility exhibits an algebraic divergence with decreasing 
temperature, replacing the logarithmic divergence found in a Fermi liquid [Yang et al, Phys. Rev. 
Lett. 106, 047004 (2011)]. In this paper we explore the single-particle and transport properties near 
the quantum critical point using high quality estimates of the self energy obtained by direct analytic 
continuation of the self energy from Continuous-Time Quantum Monte Carlo. We focus mainly on 
a van Hove singularity coming from the relatively fiat dispersion that crosses the Fermi level near 
the quantum critical filling. The fiat part of the dispersion orthogonal to the antinodal direction 
remains pinned near the Fermi level for a range of doping that increases when we include a negative 
next-near-neighbor hopping t' in the model. For comparison, we calculate the bare c?-wave pairing 
susceptibility for non-interacting models with the usual two-dimensional tight binding dispersion and 
a hypothetical quartic dispersion. We find that neither model yields a van Hove singularity that 
completely describes the critical algebraic behavior of the bare d-wave pairing susceptibility found 
in the numerical data. The resistivity, thermal conductivity, thermopower, and the Wiedemann- 
Franz Law are examined in the Fermi liquid, marginal Fermi liquid, and pseudo-gap doping regions. 
A negative next-near-neighbor hopping t' increases the doping region with marginal Fermi liquid 
character. Both T and negative t' are relevant variables for the quantum critical point, and both 
the transport and the displacement of the van Hove singularity with filling suggest that they are 
qualitatively similar in their effect. 

PACS numbers: 74.40.Kb, 71.10.Fd, 74.72.-h, 71.10.Hf 



I. INTRODUCTION 

A plausible scenario for the high temperature super- 
conductivity in cuprates is based upon the presence of a 
van Hove singularity corresponding to the saddle points 
in the single particle energy dispersionii"— These flat re- 
gions in the energy dispersion are directly observed in 
ARPES experiments on various cupratc compoundsi^"— 
Recently, it was also observed in the tunneling spectra 
of Bi-2201iiS The presence of saddle points in the energy 
dispersion is also argued to lead to a superconducting in- 
stability in other correlated systems, e.g. graphene^i. If 
the Fermi level is doped to coincide with the van Hove 
singularity, then the superconducting transition temper- 
ature can be greatly enhanced. 

The van Hove scenario is also arguedi^"— to be respon- 
sible for the marginal Fermi liquid behavior— in which 
the lifetime broadening of the quasiparticles is of the or- 
der of its energy. Thus the van Hove scenario is argued to 
account for the linear-T resistivity^i^ii^ii^ T-independent 
thermopower anomalous isotope effect?^ etc. 

There is numerical evidence for the presence of van 
Hove singularities in models of strongly correlated sys- 
tems. The energy dispersion of one hole in an antifer- 
romagnctic background has been considered in studies 
of the Hubbard mode^^^ and the t-J model These 



studies report the presence of extended saddle points. 
lAssaad and Imada ^ found that the dispersion has a 
quartic dependence with momentum near the anti-nodal 
point (tt, 0). 

These examples of extended saddle points in various 
correlated superconducting systems, and their proximity 
to the Fermi level at the doping where the maximum 
transition temperature occurs, demonstrate that it is ex- 
tremely important to understand the role played by these 
singularities. A plethora of scientific efforts have been de- 
voted towards achieving this understanding!^""— At the 
simplest level, the role of the van Hove singularity may 
be interpreted within the BCS formalism. Here, the su- 
perconducting transition temperature, Tc, is determined 
by the condition Vx'oi^ = 0) = 1, where x'o the real 
part of the q — bare pairing susceptibility, and V is 
the strength of the pairing interaction. In a BCS su- 
perconductor, Xoi^ = 0) displays a logarithmic diver- 
gence as T — )- 0, yielding the BCS exponential form for 
Tc ■ The van Hove singularity enhances the divergence 
of x'oi^ = 0), yielding higher transition temperatures. 

There is also strong evidence for a quantum critical 
point located beneath the superconducting dome in the 
cuprates, and in close proximity to the doping with the 
maximum Tc Above the quantum critical point, in 

a range of doping associated with marginal Fermi liq- 



uid behavior, the in-plane resistivity is known to vary 
hnearly with T over a wide range of teniperatures4i"— 
In the Fermi hquid region the low temperature resistiv- 
ity varies as T^. The resistivity increases as the dop- 
ing decreases from the Fermi liquid into the pseudo- 
gap region. Moreover, the thermal conductivity k,— 
the thermopower S^^^ and the tunneling conductance 
have been investigated near the quantum critical 
point of the cuprates. k, is observed to be nearly in- 
dependent of temperature in the marginal Fermi liquid 
stated and depends on 1/T in the Fermi liquid region. 
This is c onsistent with the Wiedemann-Franz Law)^ 
Kp (X T. IChakrabortv et ahl ^ suggested that the ther- 
mopower changes sign abruptly near the optimal doping 
in most of the cuprate materials, signaling a state with 
particlc-holc symmetry. Also in the marginal Fermi liq- 
uid, the tunneling conductance g{V) ^ go+gi\V\, where 
go and gi weakly depend on T and V. 

A recent study^ reported the presence of a quantum 
critical point in the two-dimensional Hubbard model, 
where the quasiparticle spectral weight becomes zero. 
This quantum critical point separates the non-Fermi liq- 
uid pseudogap from the Fermi liquid region, and is sur- 
rounded by a superconducting dome (c.f. the inset in 
Fig. [T|). At finite temperatures, the Fermi liquid and 
pseudogap regions are separated by the marginal Fermi 
liquid. Interestingly, at the quantum critical point, the 
density of states is found to be nearly particle-hole sym- 
metric at low frequencies with a sharp peak at = 0. 
This filling is tantalizingly close to the optimal dop- 
ing where the superconducting transition temperature 
Tc attains its maximum. The proximity of the super- 
conducting dome to th e quantum c ritical point was re- 
cently investigated by lYang et alj ^. Unlike the BCS 
case, they found that the bare d-wave pairing suscepti- 
bility x'odi'^ — 0) diverges algebraically as — = at the 

quantum critical point, thus leading to a strongly en- 
hanced Tc- Using the Kramers-Kronig relation between 
the real part and the imaginary part of the susceptibil- 
ity, Xod(^) — ~ y ^^'''^ — ^dw, the algebraic divergence 
of XodC^) found to come from a scaling behavior of 
the imaginary part x'odi^)- When T^^'^x'odi^)/^ plot- 
ted against w/T, the different temperature curves fall 
on top of each other when uj > Ts = AtT/J, determin- 
ing a scaling function H{x) such that T'^^^x'odi^)/'^ ^ 
H{uj/T) w (w/T)-3/2 (see Fig. [T]). The contribution 

y-3/2 r 

from H to x'{T) = —— / H{uj/T)duj cx T^^/^ which 

will dominate at low T . 

Since this enhanced behavior is expressed in the bare 
pairing bubble, dressed by the self energy but with no 
vertex corrections, this discussion naturally raises the 
question about the role played by the van Hove singu- 
larity in the quantum criticality and its possible connec- 
tion to the superconducting T^. In this manuscript, we 
use the dynamical cluster quantum Monte Carlo method 
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FIG. 1: (color online) Frequency dependence of the imag- 
inary part of the particle-particle d-wave susceptibility ob- 
tained from the Dynamical Cluster Approximation for vari- 
ous temperatures at the critical filling n = 0.85 when U = 6t, 
t' = Q and At = 1. The arrow denotes the direction of 
decreasing temperature. All the curves fall on top of each 
other for frequencies greater than Ts/T ~ 4f/J. The inset 
shows the phase diagram for the same parameters {t, t' and 
U) including the Fermi Liquid (FL), Marginal Fermi Liquid 
(MFL) and Pseudogap (PG) regions. The lines indicate the d- 
wave transition temperature Tc (determined by extrapolation 
of data from A^c = 8, 12 and 16 site clusters), the pseudo- 
gap tempera ture T* , and the Fermi liquid temperature Tf- 
(Taken from lYang et al.P ^) The quasiparticle fraction on the 
Fermi surface vanishes at the quantum critical point where 
T* — Tf — 0, and remains zero in the pseudogap region.— 



to explore the relationship between the quantum critical 
point and the van Hove singularity for high-temperature 
superconductivity in the Hubbard model. We obtain high 
quality estimates of the real- frequency single-particle self 
energy at the cluster momenta K by direct an- 

alytic continuation of the Matsubara frequency self en- 
ergy S(iir, iujn) using the maximum entropy methodi^i^ 
This direct method avoids the artifacts on the self en- 
ergy that come about by inverting the coarse-grained 
Dyson's equation]^ In the model without next-near- 
neighbor hopping {t' = 0), we find that, as we dope 
the system across the quantum critical filling, a flat re- 
gion in the dispersion crosses the Fermi level, accompa- 
nied by a sharp nearly symmetric peak in the density of 
states which also passes through the Fermi level. We find 
that the resistivity follows a linear-T dependence over a 
wide range of temperatures yet a narrow range of dop- 
ing (sec Fig. [T4|) . We use these high quality estimates 
of the self energy to calculate the bare pairing suscep- 
tibility, we again find the collapse of the data found in 
Fig. [TJ To understand the role played by the van Hove 
singularity in determining this critical behavior, we have 
calculated the pairing susceptibility in the d-channel for 
two non-interacting models at half filling - the standard 
quadratic dispersion and a hypothetical quartic disper- 
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sion. While the quartic dispersion can yield the observed 
algebraic divergence of Xod('^ = 0), neither dispersion 
produces the collapse of the data found in Fig.[Tl suggest- 
ing that a van Hove singularity alone does not capture 
this phenomena. For negative t', the resistivity follows 
a linear-T behavior over a wider range of doping, and 
the sharp peak in the density of states and the flat re- 
gion of the dispersion linger near the Fermi level for the 
same wider range of doping. These results suggest that 
the doping region affected by quantum criticality at low 
temperature becomes larger when t' < 0. We also show 
that the zero-frequency imaginary part of the self energy 
Yi"{T,uj = 0), the dominant contribution to the resistiv- 
ity, has a wider range of linear-T behavior for t' < than 
t' = Q. All this motivates us to speculate a phase dia- 
gram near the quantum critical point in the discussion 
section. 

This paper has been organized as follows. Section [IT] 
briefly outlines the model and methods used in this study. 
Results are presented in section Hill Single particle prop- 
erties are discussed in section HlIAl the pairing suscepti- 
bility calculation in section HlI Bl the effect of t' on the 
dispersion in IIII C[ and transport results in IIII Dl The 
broader implications of our results are discussed in sec- 
tion [IV] and the paper is concluded in section |Vl 

II. FORMALISM 

In this work, we look for direct evidence of the van 
Hove singularity and marginal Fermi liquid behavior in 
the spectra, electronic dispersion, and transport proper- 
ties of the two-dimensional Hubbard model 

ka ' 

where ^^ka^ ^^"^ creation (annihilation) operator 

for electrons with wavevector k and spin cr, rijo- = cl^Cia 
is the number operator, and the bare dispersion is given 
by 

= — 2t (cos kx + cos ky) — At' (cos k^ cos ky — 1) , (2) 

with t and <' being the hopping amplitude between the 
nearest-neighbor and the next-near-neighbor sites, re- 
spectively, and U is the on-site Coulomb repulsion. 

We employ the Dynamical Cluster Approximation 
(DCA)^iS^ with a quantum Monte Carlo algorithm as 
the cluster solver. The DCA is a cluster mean-field theory 
that maps the original lattice onto a periodic cluster of 
size Nc = embedded in a self-consistently determined 
host. This many-to-one map is accomplished by dividing 
the lattice Brillouin zone into cells centered at momenta 
K, and coarse graining the lattice Green's functions by 
summing over the momenta labeled with k within each 
cell 

G(fs:,c.) = ^^G(ii: + k,c.), (3) 

k 



where G and G are the coarse-grained and the lattice 
single-particle propagators, respectively. The coarse- 
grained Green's function defines the cluster problem. 
Spatial correlations up to a range Lc within the clus- 
ter are treated explicitly, while those at longer length 
scales are described at the mean-field level. However 
the correlations in time, essential for quantum criticality, 
are treated explicitly for all cluster sizes. To solve the 
cluster problem, we use continuous-time quantum Monte 
Carlo^, which has no Trotter error,— and the Hirsch-Fye 
quantum Monte Carlo method^i^ for the charge polar- 
izability in Fig. [5] We employ the maximum entropy 
method^ to calculate the real-frequency spectra. 



A. Calculation of Single-Particle Spectra 

In previous calculations of the single-particle spec- 
tra, we analytically continue the quantum Monte Carlo 
G{K,t) to obtain G{K,uj), and then invert the coarse- 
graining Eq. ^ to obtain the self energy T,{K,co). This 
last step can introduce spurious features in E(if,u;). As 
observed previously?^ it is better to analytically continue 
the self energy directly. However, the self energy spectra 
does not share the normalization of / dijjA{K ,uj) = 1, 

where A{K,lj) = G" {K This normalization is a 

TT 

desirable feature since it allows us to treat the spectrum 
as a normalized probability distribution. Since the Hub- 
bard self energy T,{K,iuJn) = T,h + U'^X<y,ali^n + • ■ ■ , 
where x<T,cr = {nana) — {na)"^ = '^(t(1 — n„) is the local 
polarizability of a single spin species cr, and S(i<r, iLOn) — 

T,H = / duj — ^ — ^ — ^. It is easy to see that the in- 

TT J iuj„ — iO 

tegral of J^{K,iu!n) — is U'^Xa,a- Therefore we will 
analytically continue 

i:{K,iuJn) -JIh f, <j{K,uj) 

— ^ / duj- , (4) 

U''Xa,a J UVn-U! 

where a{K,uj) = --j:"{K,uj)/U^Xa.a, Jduja{K,uj) = 

TT 

1, using Xa.a calculated in the Monte Carlo process. Af- 
ter that we obtain the lattice self energy J^{k,uj) by in- 
terpolating the cluster self energy Y,{K,uj) to get the 
single-particle spectral function A{k,uj). 



B. d-wave Pairing Susceptibility 

We calculate the susceptibility in the c?-wave channel 
to the pair field V = ~fdi>\ + h.c, for various models 
with a van Hove singularity at the Fermi level. Here 

^rf = 5 Si {^\+x ^ singlet creation operator, 

where h\^^ creates a singlet at bond -t- d), a ^ x,y, 
and fd is a complex constant. The non-interacting d-wave 
pairing susceptibility xod can be computed by calculating 
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the susceptibility bubble 

Xod(T) ^tJ2 5d(fe)G°(fc, *L.„)G'°(-fe, -jc^„), (5) 

where gd{k) is the d-wave form factor given by gd{k) = 
cos fcj; — cosfcj,. G'~'{k,iu!n) is the non-interacting Green 
function given by 

G"(fc,»^«)- . ^ ■ (6) 

with the bare band dispersion in Eq. ([2]). xod can be 
evaluated using standard Matsubara summation which 
gives 

Xo4T)=E5d(fc)(^^|^y (7) 
where is the Fermi function. 

C. Transport Coefficients 




CO 



FIG. 2: (Color online) The local r = imaginary part of the 
dynamical charge polarizability divided by the initial slope at 
= 0, for n = 0.85, U = 6t, At = 1, t' = and Nc = 16. 
It satisfies the marginal Fermi liquid form given by Eq. 
Inset: the zero frequency slope of x"(^ = 0,cj) is roughly 
linear in inverse temperature, as expected. The line is a linear 
fit to the expression a + b/T. 



To explain the anomalous transport properties of the 
marginal Fermi liquid, Varma et aL— postulate that 
for a wide range of wavevectors q, excitations make a 
contribution to the absorptive spin and charge polariz- 
abilities reflected by 



X"iq, w) cx min {\uj/T\ , 1) sign(w). 



(8) 



Electrons scattering from these excitations acquire a self 
energy 



S(fc, uj) oc Lu In (x/iOc) — i'Kx/2, 



(9) 



where x = max(|ti;|, T), and Wc is a cutoff. The marginal 
Fermi liquid ansatz has several consequences on experi- 
mentally relevant quantities, including transport anoma- 
lies, such as the linear-T electrical resistivity, the tun- 
neling conductance giV) ~ go + gi\V\, the photoemis- 
sion, the nuclear relaxation rate T^^ ^ aT + b, the op- 
tical conductivity <y{i-o), the Raman scattering, and the 
superconductive pairing. For the specific heat Cy{T) 
and thermal conductivity k(T), Varma argued that the 
normal state's electronic contribution is hard to extract 
from the experimental data due to the large phonon con- 
tribution. The electronic thermal conductivity for the 
marginal Fermi liquid approximates to a constant be- 
cause the Wiedemann-Franz law roughly holds. 

To calculate the various Onsagcr transport coefficients 
we use the Kubo formula)^ 



where v°'[k) is the a-componcnt of the electron group ve- 
locity and A(k, uj) is the single-particle spectral function. 
The different transport coefficients are given by combina- 
tions of . For example, in units where e = 1 and the 
chemical potential = 0, the resistivity p{T) = 
the thermopower S = —L^'^/TL^^, the thermal conduc- 
tivity 7^ {L^^ - {L^^f/L^^), and the Pcfiier coeffi- 
cient n = L^i/iii. 

We note that a simpler estimate exists for the ther- 
mopower S. Here, we perform a Sommerfeld expansion 
of L^^ at the Fermi level and get an alternative form: 



3 duj 



LJ = 0' 



(12) 



If the electron group velocity is a constant, and the 
square of the single-particle spectra is approximated by 
5{u) — ^k)^k^ where tj^ is the relaxation time, also as- 
sumed constant, the thermopower over temperature be- 
comes just the derivative of the logarithm of the density 
of states at the Fermi levelj2£ 



III. RESULTS 



A. Single Particle Properties for t' = Q 



L% = 7tJ dco oj^+^-' I?„^(c.), (10) 

where / is the Fermi function and 

^"/5H = ^E^"(^>^W^(fc'^)'' (11) 

k 



We first explore the charge polarizability in the 
marginal Fermi liquid region at rt = 0.85. The imaginary 
component of the local charge polarizability Xc{r ~ 0, a;) 
is plotted in Fig [2j The main plot shows x"(^ = 0,uj) 
divided by its initial slope at zero frequency (determined 
in the inset), so that the curves coincide for low w. 
At higher frequencies, the curves break from this linear 
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FIG. 3: (Color online) Frequency dependence of the imagi- 
nary part of the self energy at the Fermi momenta, E"(A;f, i^), 
along FM and XF for U = 6t, At = 1, t' = 0, Nc = 16 and 
13 — 58. Right panels show a zoom of the low frequency re- 
gion. Dashed lines fit the data linearly for n = 0.85 and 
quadratically for n = 0.75. 
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FIG. 4: (Color online) Temperature dependence of the imagi- 
nary part of the self energy T,"{kF,ijJ = 0) at the Fermi energy 
and momenta, U = 6t, 4:t = 1, t' — and Nc — 16. The linear 
dashed lines fit the self energies for n — 0.85 and 0.86 below 
T = 0.031. 



rise at a frequency roughly proportional to the temper- 
ature. The inset shows that the zero frequency slope, 
Xcif = Oi'^)/<^L=o' roughly linear in inverse temper- 
ature up to T « 0.2 or roughly 2 J = 8t'^/U. These 
features are consistent with the marginal Fermi liquid po- 
larizability in Eq. ([8]) . The spin polarization (not shown) 
docs not display such an extended region of marginal 
Fermi liquid character. This result is consistent with 
marginal Fermi liquid behavior due to the proximity of a 
quantum critical point associated with phase separation. 

Fig. [3] shows the frequency dependence of the imagi- 
nary self energy at the Fermi momenta along the anti- 



FIG. 5: (Color online) The single-particle density of states for 
t' ^ 0, U = 6t, At = 1, Nc = 16 and /3 = 58. The density of 
states shows low-energy particle-hole symmetry at the filling 
of n = 0.88. 



nodal (Xr) and the nodal (FM) directions for three fill- 
ings: n = 0.75 (Fermi liquid), n = 0.85 (marginal Fermi 
liquid) and n = 0.95 (pseudo-gap). For the self energy at 
n = 0.75, the quadratic dashed line (bottom right panel) 
provides a good fit to S"(fei?,a;) for small w, as expected 
from the Fermi liquid theory.— The marginal Fermi liq- 
uid self energy has a form given by Eq. ^ , which states 
that the imaginary self energy is proportional to the neg- 
ative temperature for small frequency and to the nega- 
tive uj when the temperature is low. The marginal Fermi 
liquid self energy in Fig. [3] (upper right panel) shows a 
linear behavior, but interestingly with different slopes 
for positive and negative w. This is not consistent with 
Eqn. ([3]), but this may b e due to the pres ence of some 
short range order. I.e., iMarkiewicz et all ^ calculated 
the self-energy due to the random-phase approximation 
(RPA) magnetic polarizability for a single band model 
with dispersion obtained by a fitting to the tight bind- 
ing local-density approximation. They found that the 
self-energy has linear forms but with different slopes on 
positive and negative w when the van Hove singularity is 
at the Fermi level and quadratic otherwise. 

Fig. m shows the temperature dependence of the self 
energy when w = and t' = 0. Again we find a result 
consistent with Eqn. ([S]), I]"(fcj7,a; — 0) oc — T, around 
the marginal Fermi liquid filling for n = 0.85 and 0.86. 
The dashed lines are linear fits. The error bars are esti- 
mated by changing the random seeds in the Monte Carlo 
calculation. However, they do not refiect the system- 
atic error that comes from the bias towards the default 
model, which in this case is the spectra from the next 
higher temperature. This error is largest at low T where 
the data is weak due to the minus sign problem. 

The self energy and dynamic charge polarizability near 
the critical filling Uc = 0.85 arc generally consistent with 
Varma's marginal Fermi liquid ansatz, as are the results 
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found previously for the kinetic and potential energies 
which vary with temperature like ln(T)S and the van- 
ishing wave function renormalization factor—. To under- 
stand the relationship of these results to the van Hove 
singularity, we will explore the density of states and the 
quasiparticle dispersion. 

The density of states for several fillings is shown in 
Fig. [5l Since we have highly enhanced the quality of the 
self energy by direct analytic continuation of E(/C, 
the density of states in Fig. [3 shows sharper featu res as 
compared to the results of Vidhvadhiraia et alj ^. As 
the doping decreases from the Fermi liquid to the pseudo- 
gap region, the peak of the density of states moves from 
positive to negative energy while its intensity is reduced. 
For n — 0.95, a pseudogap begins to open and a peak to 
form at positive frequencies. The half- filled case (n = 1, 
not shown) shows upper and lower Hubbard bands lo- 
cated at positive and negative frequencies, respectively. 
The density of states for filling 0.88, close to the critical 
filling of 0.85, shows low-frequency particle-hole symme- 
try. 

Fig. [6] shows the dispersion obtained from the peaks 
of the spectral function for four fillings: n = 

0.85, 0.87, 0.88 and 0.95, along the anti-nodal direc- 
tion and around the Fermi vector kp. In order to de- 
fine k p , we look for the maximum value of the zero fre- 
quency spectral function A{k^uj = 0). In the Fermi liq- 
uid and marginal Fermi liquid regions, this definition is 
roughly equivalent to the Luttinger surface defined where 
G'{k,uj = 0) changes sign. The two definitions yield dif- 
ferent results for the pseudogap region, especially when 
t' < which enhances the pseudogap. However, this dif- 
ference is not large enough to qualitatively change our 
results or to change any of our conclusions. Therefore, 
for simplicity, we only show results using the first defini- 
tion of kp. For a particular filling, the left panel shows 
the dispersion along the direction [ky — 0), while the 
right panel shows the dispersion along ky (k^ = kp^). 
A common identifiable feature for all fillings is the pres- 
ence of a flat region in the dispersion. This flat region 
is responsible for the van Hove singularity in the density 
of states. The van Hove singularity passes through the 
Fermi level at a filling of ri ss 0.88, which is near the quan- 
tum critical filling, n sa 0.85, where the quasi-particle 
weight Z goes to zero.— At this filling, the topology of the 
Fermi surface also changes from hole-like (closed around 
the point k = {tt,tt)) to electron-like (closed around 
k = (0,0)) with increasing filling (not shown) as seen 
in experiments.— The dispersion along the ky direction 
remains pinned near the Fermi level for a range of dop- 
ing near the center of the superconducting dome, while 
the dispersion along the k^ direction passes continuously 
through the Fermi level. This anisotropic motion of the 
flat dispersion is consistent with a van Hove peak which 
moves continuously through the Fermi level as shown in 
Fig. [21 and would correspond to a flat region at the Fermi 
level which is most isotropic at the crossing and shrinks 
to narrow pencil-like regions for fillings above and below 



the crossing. 

The dispersion along the anti-nodal direction as a func- 
tion of kx for various fillings is displayed in Fig. [71 In- 
terestingly, a quadratic form fits well to the data for all 
fillings. Next, we investigate if such a dispersion can 
capture the critical algebraic divergence of the pairing 
susceptibility. 



B. Pairing Susceptibility 

The density of states and the dispersion show clear evi- 
dence for a van Hove singularity which crosses the Fermi 
level near the critical filling. In order to see whether 
the van Hove singularity alone is sufficient to explain the 
enhanced bare pairing bubble, we calculate the pairing 
susceptibility in the d-wave channel for two simple mod- 
els having a van Hove singularity at the Fermi level. We 
begin with the tight binding model given by Eq. ([5]) at 
half filling and t' = 0. The associated density of states 
has a logarithmic singularity at e = 0, N{e) = log|e|. 
The temperature dependence of Xod '^^^ be obtained by 
converting Eq. ([7]) to an integral over energy with a tem- 
perature T cutoff. It results in a — (logT)^ behavior. 
This is confirmed by explicit calculation of the sum in 
Eq. ([7]) as illustrated in Fig. [H As shown in the inset, 
the imaginary part does not show scaling behavior as 
seen in Fig. [H 

We also consider the next higher order model allowed 
by the symmetry of the square lattice, a hypothetical 
model with a quartic dispersion 



= --((|fc.|-^)^-fc^4). 



(13) 



Such an extended form has been observed in experi- 
ments^ and also confirmed by theoretical studiesJ^ The 
low energy density of states for the quartic dispersion 
becomes N{e) ~ l/y^.^l Following a similar logic as 
we used for the tight binding dispersion, for a quartic 

dispersion, we get x'od ~ 

Results for the explicit calculation (Eq. ([7])) are shown 
in Fig. [SI and are consistent with the analytical argu- 
ments above. Though the temperature dependence of 
the real part of the bare susceptibility is found to be al- 
gebraic for this quartic dispersion, the inset reveals that 
the imaginary pai ring suscept ibility does not exhibit the 
scaling found by lYang et al.l ^. Thus, the simple non- 
interacting picture of the van Hove singularity at the 
Fermi level does not completely describe the true tem- 
perature and frequency dependence of the susceptibility 



C. Effect of negative t' 

The single-band Hamiltonian used to model the hole- 
doped cupratcs generally includes a negative next-near- 
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FIG. 6: (color online) Energy dispersion obtained from the peaks of the spectral function A{k, oj) for various fillings around 
the Fermi vector kp along the anti-nodal direction for t' = 0, U = 6t, At = 1, Nc = 16 and /3 = 58. By fixing fcj, = we explore 
the dispersion along the kx direction, and for kx = kpx the dispersion along the ky direction is plotted. Notice that the energy 
axes are inverted so that positive energies are plotted down. The dispersion along the ky direction remains pinned near the 
Fermi level for a range of doping near the center of the superconducting dome (c.f. the inset of Fig. [l| 

neighbor hopping t' . For t' = — O.lt, the temperature 
dependence of the self energy at the Fermi momenta and 
energy is shown in Fig. (TU] We find that J^"{kF,LL! — 0) 
follows a linear behavior over a wider range of fillings, 
from n = 0.83 to 0.87. 

Fig. [11] demonstrates that the inclusion of a negative 
t' also results in the pinning of the flat part of the ky- 
dispersion to the Fermi level. However, now the pinning 
is observed for a larger range of fillings, roughly 0.80 
to 0.86. Thus both measurements, the temperature de- 
pendence of the self energy and the pinning of the flat 
dispersion to the Fermi level, are consistent. If we take 
the viewpoint that the quantum critical point and the 
pinning of the dispersion along ky to the Fermi level are 
concomitant aspects of quantum criticality, then a nega- 
tive t' leads to a larger range of quantum critical fillings. 
We will also see the signature of this behavior in various 
transport properties discussed in Section UlI Dl 

Fig. [T2l shows the density of states for t' /t = —0.1 and 
various fillings as a function of w. For a given filling, the 
inset of FigfT^shows that the peak in the density of states 
is slightly shifted to smaller frequencies when compared 
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FIG. 7: (color online) Single particle dispersion around the 
Fermi energy taken along the anti-nodal direction {ky = 0). 
The data is from Fig. [6l The solid lines are fits to a quadratic 
dispersion. 
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FIG. 8: (color online) Temperature dependence of the real 
part of the particle-particle d-wave susceptibility at cj = 
for the two-dimensional tight binding dispersion at half fill- 
ing. Note that Xod diverges logarithmically as T — >■ 0. Inset: 
Frequency dependence of the imaginary part of the particle- 
particle d-wave susceptibility. Note that the curves corre- 
sponding to various temperatures do not scale at large fre- 
quency. The arrow denotes the direction of decreasing tem- 
perature. 
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FIG. 9: (color online) Temperature dependence of the real 
part of the particle-particle d-wave susceptibility ak, uj = Q for 
the quartic dispersion of Eq. (|13|) at half filling. Note that x'od 
diverges algebraically ~ as T — ^ 0. A fit to a + h/x"^ 

gives values of a = —10.6, b = 1.98 and c = 0.54. Inset; 
Frequency dependence of the imaginary part of the pairing 
d-wave susceptibility. Note that the curves corresponding to 
various temperatures does not scale well at large frequency. 
The arrow denotes the direction of decreasing temperature. 



with the peak in the density of states for i' = 0. It 
displays particle-hole symmetry roughly at n = 0.84, not 
n = 0.88 as for t' = 0. Moreover, if we use Awp/An, 
where Wp is the location of peak in the density of states 
and n is the fihing, to estimate the rate at which the 
peak crosses the Fermi level, we find that the peak of the 
density of states for t' = crosses the Fermi level more 




FIG. 10: (Color online) Temperature dependence of the imag- 
inary part of the self energy E"(fcF,i^ = 0) at the Fermi en- 
ergy and momenta, U = 6t, Nc = 16, 4t = 1, and t'/t = —0.1. 
The self energy for filling between n = 0.83 and 0.87 shows a 
linear-T behavior. 



quickly than the peak for negative t' . This can be seen 
in the inset of Fig. [121 where the filling dependence of 
the peak location has a steeper slope for t' = than that 
for t'/t = —0.1 at the Fermi level. This confirms that 
negative t' leads to a wider range of fillings with a van 
Hove peak near the Fermi level. The fact that this range 
of fillings coincides with the region where marginal Fermi 
liquid behavior is seen in the self energy suggests that the 
van Hove singularity and quantum criticality are related. 

Another interesting point to be noted here is that, 
when compared to the t' = result, the quasi-particle 
peaks become more incoherent for negative t' . This 
can be seen in the Matsubara quasiparticle weight along 
the antinodal momentum direction, Zan^ displayed in 
Fig. [13] as a function of temperature for different fill- 
ings. The c^uasiparticlc fraction is consistently smaller 
for t'/t = -0.1 than for t' = for all fillings. This can 
also be seen through the increase of the blue color in the 
dispersion curves in Fig. [11] when compared with Fig. [6] 



D. Transport Properties 

Matrix element effects^^i^, and the low precision of in- 
verse photoemission can complicate the direct measure- 
ment of the flat dispersion resulting in the van Hove sin- 
gularity, making indirect probes like the Fermi surface 
topoIogy2iiZl and transport measurements more impor- 
tant. The van Hove singularity and the quantum critical 
point will also impact the transport properties of the sys- 
tem. Using the Kubo formula under the relaxation time 
approximation in Eqs. (|10[) and (|11[) . we obtain the re- 
sistivity, thermal conductivity, and thermopower in the 
Fermi liquid, marginal Fermi liquid, and pseudo-gap re- 
gions. 
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FIG. 11: (color online) Energy dispersion obtained from the peaks of the spectral function A{k,ijj) for various fillings around 
the Fermi vector fe^ along the anti-nodal direction for t' /t = —0.1, U = 6t, 4t = 1, Nc — 16 and /3 = 48. The dispersion 
along the ky direction remains pinned near the Fermi level for a broader range of dopings than found when t' — 0. Again, the 
dispersion along the kx direction moves continuously across the Fermi level. 
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FIG. 12: (Color online) Single-particle density of states for 
t'/t = -0.1 when U = 6t, = 16, 4t = 1, and /? = 48. Inset: 
Comparison of the filling dependence of the position of the 
peak of the density of states {cjp) for t' = and t'/t = —0.1. 
As the filling changes, ujp for t' = crosses the Fermi level 
more quickly than for t'/t — —0.1. 




FIG. 13: (Color online) Temperature dependence of the Mat- 
subara fraction along the anti- nodal direction, Zan, for vari- 
ous fillings for t' — and t'/t — —0.1. At the same filling and 
temperature, the Matsubara fraction decreases when t' < 0. 



Fig. [T4l shows the resistivity as a function of tempera- 
ture for t' — 0, left panel, and t'/t = —0.1, right panel. 
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FIG. 14: (Color online) Resistivity versus temperature for t' — (left) and t' /t = —0.1 (right) with U = 6t, Nc = 16, and 
At = 1 . The dashed lines in the insets are linear fits. For t' = 0, the resistivity shows linear-T behavior for n = 0.85 and 
n = 0.86. For t' /t — —0.1, the resistivity shows linear-T behavior from n = 0.83 to n = 0.87. 



Linear resistivity reveals evidence of the marginal Fermi 
liquid because the electronic cross section is proportional 
to —T,"{kp,uj ~ 0) at low T, and, as seen in Fig. Uand 
[TOl this self energy is linear in T. Again, for t' = a 
narrovs^ range of fillings, from n = 0.85 to 0.86, displays 
a linear-T resistivity at low T. While for t' /t = —0.1 a 
larger range of filling, n = 0.83 to 0.87, exhibits a linear 
temperature dependence. The linear resistivity in the 
marginal Fermi liquid region is consistent with experi- 
ments.ii."il For n = 0.75, both t' = Q and t' /t = -0.1 
show Fermi liquid character, with a resistivity which goes 
to zero quadratically when T approaches zero. The fact 
that doping region with marginal Fermi liquid character 
increases with negative t' has consequences for the phase 
diagram near the quantum critical point, which we will 
discuss in Section ITVl 

According to the Wiedemann- Franz Law, k/{(jT) = 
7r^/3 (ks = e = 1), the thermal conductivity of a Fermi 
liquid is inversely proportional to Fig. [T5] shows 

that K (X l/T for n = 0.75 when t' /t = and -0.1, , 
but weakly depends on T for the marginal Fermi liquid 
and the pseudo-gap regions. The inset shows that, for 
n = 0.75, the Wiedemann- Franz ratio n/ (aT) approaches 
a constant which is less than tt^/S when T < 0.08. 
iDahm et aH ^ investigated the two-dimensional Hubbard 
model for n ~ 0.9 and also found a smaller Wiedemann- 
Franz ratio. However, we find that the Wiedcmann-Franz 
ratio is larger than 7r^/3 for the marginal Fermi liquid 
{n = 0.85) and pseudo-gap (n ~ 0.95) regions. We also 
see that the thermal conductivity becomes very small 
as T — > for n = 0.95 and saturates to a constant for 
n ~ 0.85. So, when studying k, the marginal Fermi liquid 
seems to separate the Fermi liquid from the pseudo-gap 
region. The dashed curves in Fig.[Tn]for t' /t — —0.1 data 
are always below the solid curves for t' — Q when plot- 
ting K. However, the t' jt = —0.1 data is above the = 
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FIG. 15: (Color online) Thermal conductivity versus temper- 
ature for U = 6t, Nc = 16, and At = 1. Inset: Wiedemann 
Franz ratio for the same physical parameters. The horizontal 
solid line labels the constant tt'^/S. 

results when we focus on the ratio k/uT. This implies 
that negative t' reduces the electrical conductivity more 
than the thermal conductivity. 

IChakrabortv et alT ^ argue that the thermopower 
changes sign near the quantum critical point, and that 
this is related with the development of a state with 
particle-hole symmetry. Fig. [12] shows the thermopower 
S as a function of filling. For t' = and /3 = 58, the filling 
at which 5* changes sign is roughly 0.80. We expect that 
the zero-crossing of the thermopower will approach the 
critical filling of 0.85 for decreasing T. However, this is 
different from the filhng, n = 0.88, at which the density 
of states displays particle-hole symmetry. 

Wc find that due to the fc-depcndcncc of the relax- 
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FIG. 16: (Color online) Thermopower as a function of filling 
for U = 6t, Nc = 16, At — 1. Lines are guides to the eyes. 
Inset: Frequency dependence of ©^(a;) (c.f. Eq. (jlip for 
different fillings when t' = 0, p = 58. The slope of D^xiij^) 
at a; = is proportional to the thermopower according to 
Eq. dni). 



ation time and the electron group velocity, the filling at 
which the thermopower crosses zero does not occur at 
the filling where the density of states shows a particle- 
hole symmetry at low energies. In Fermi liquid theory^ 
if we assume constant relaxation time and group velocity, 
the thermopower is proportional to the derivative of the 
logarithm of the density of states at the Fermi level. This 
would suggest a thermopower which changes sign as the 
van Hove singularity crosses the Fermi level. However, 
in this approach ^(fc, w)^ in Eq. (|lip is approximated by 
S{lli — (■if)T, where r = is a fc-independent relaxation 
time. In addition, the electron group velocity we use in 
our calculation also has a fc-dcpendence: 

w^(fc) = = 2tsmkx + W sin cos ky. (14) 

If we compare the quantity V^xi^) (the inset of Fig fTB)) 
and the density of states (Fig. [S]) for different fillings and 
the same t' = 0, we find that the effect otv^{k)^ is to pull 
the peak of the density of states to the left, because the 
sin kx term suppresses the contribution of the van Hove 
singularity at X{tt, 0) and enhances the contribution from 
the states below the Fermi level. As a result, the ther- 
mopower changes sign continuously and at a filling where 
the quantity I?(a;) has a zero slope at the Ferm i level, 
around n = 0.85. As noted by IChakrabortv et ahl ^. this 
filling is different from the one where the density of states 
displays particle-hole symmetry, n = 0.88 for t' = 0. The 
impact of the van Hove singularity on the thermopower 
and other transport coefficients is diminished by the fact 
that the van Hove singularity comes predominantly from 
a region in k space where the group velocity goes through 
zero. 



IV. DISCUSSION 

The results presented here have implications for both 
the quantum critical phase diagram and the proximity of 
the superconducting dome to the quantum critical point. 

The close proximity of the quantum critical filling and 
the filling where the van Hove singularity crosses the 
Fermi level suggests that the two are related or even con- 
comitant. Near the quantum critical point, we find that 
the fiat part of the dispersion orthogonal to the antinodal 
direction is pinned to the Fermi level, but not the disper- 
sion along that direction (Fig. [5]). We also find that the 
low energy density of states exhibits particle-hole sym- 
metry (Fig. [5]). The linear-T resistivity and self energy, 
characteristic of a marginal Fermi liquid, are observed for 
the same fillings where this pinning is observed. Within 
the Dynamical Mean Field Approximation (DMFA)i^, it 
is known that if the van Hove singularity is pinned to 
the Fermi level for the non-interacting case, it remains 
pinned even for the interacting case due to the momen- 
tum independent nature of the self-energy. In addition, a 
van Hove singularity initially away from the Fermi level 
will tend to move towards the Fermi level due to the nar- 
rowing of the coherent component of the band resulting 
from electronic correlations. In the simplest Fermi liquid 
picture, the coherent part of the single-particle Green 
function G{k,uj) = Z{k)/{uj + iQ+ - Z{k)e{k)). So, ff 
Z{k) becomes small for some values of fc, then we would 
expect to see a flattening of the observed quasiparticle 
dispersion, Z{k)e{k), accompanied by a shift of the peak 
towards the Fermi level with a vanishing weight. The new 
result of our work is that the non-local correlations in- 
cluded in the DCA, but absent in the DMFA, are able to 
move the van Hove singularity, with finite weight, to and 
even across the Fermi level. This cannot be due solely 
to the narrowing of the coherent part of the band, since 
the van Hove singularity crosses the Fermi level where 
quasiparticle fraction Z is already zero. 

Neither a non-interacting picture of the van Hove sin- 
gularity can completely describe the superconducting 
transition in the vicinity of the critical filling. In a BCS 
superconductor, the transition is driven by a logarithmic 
divergence of the bare pairing bubble as the temperature 
falls. In a recent work^ the bare d-wave pairing suscep- 
tibility x'od of the 2D Hubbard model was found to di- 
verge algebraically as — ^ at the quantum critical filling, 

instead of logarithmically, giving rise to a higher Tc . In 
the simulation, we traced the origin of this algebraic be- 
havior to a component of the dynamic bare bubble which 
scaled as Xodi^)/^^ = T-^^'^H{uj/T) (see Fig.[T]). A van 
Hove singularity is known to enhance the divergence of 
the bare pairing bubble. Since the bare d-wave pairing 
susceptibility is dressed only by the self energy, with no 
vertex corrections, a van Hove singularity seems to be the 
most likely explanation of its enhanced divergence. How- 
ever, we found that a simple non-interacting picture with 
a van Hove singularity at the Fermi level does not com- 
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FIG. 17: Schematic chemical potential-temperature (/i — T) 
phase diagram for three values of t': (from left to right) t' > 0, 
t' — 0, t' < scenario presented in this paper and a t' < 
scenario based on a Lifshitz transition. CC (QC) indicates 
the classical (quantum) critical region. ML (MG) indicates 
the Mott liquid (gas) region. T* is the crossover temperature 
between the Mott liquid and the critical region, Tfl sepa- 
rates the Mott gas from the critical region. Note here, we 
have ignored other phases to focus attention on the quantum 
critical region. 



pletely explain the observed phenomena. The standard 
quadratic dispersion gives a logarithmic divergence of x'od 
for the half filled model. A hypothetical quartic disper- 
sion yields the observed algebraic divergence for Xod^ 
does not give the correct scaling for the imagin ary part 
of the bare susceptibility found in lYang et alj ^. The 
latter is a consequence of the proximity to a quantum 
critical point, but not necessarily part of the van Hove 
scenario. 

Our results also shed additional light on the quan- 
tum critical phase diagram. We found previousl y^^i^" 
that when t' > 0, there is a first order phase separa- 
tion transition with a first-order line of co-existence in 
the /i — r phase diagram as shown in Fig. 1171 In anal- 
ogy with liquid-gas mixtures, we identify the two phases 
as a Mott liquid, which is insulating and incompressible 
with well formed local moments and short ranged order, 
and a Mott gas, which is a weakly compressible metallic 
Fermi liquid. The two phases have the same symme- 
try, so the first order line of co-existence terminates at 
a second order point where the charge susceptibility di- 
verges. In analogy to other liquid-gas mixtures, there is a 
fan-shaped region dominated by fluctuations for temper- 
atures above the critical point with neither liquid nor gas 
character. When t' — 0, the critical point goes continu- 
ously to zero temperature and thus becomes a quantum 
critical point2£. Above the quantum critical point, the 
marginal Fermi liquid region is found to exist in the V- 
shape quantum critical region. Inside this region, the 
only scale is the temperature (like Eq. ^ and T* 
is the temperature separating the marginal Fermi liquid 
from the pseudo-gap region. T* does not separate the 
quantum critical region from a region of hidden order. 
Rather in our scenario, it is only the boundary of the 
quantum critical region. As we cross from the quantum 
critical region to the Mott liquid, the character of the 
Mott liquid becomes apparent including a pseudo-gap. 

Here we consider the effect of a negative next-near- 
neighbor hopping t' on the single particle dispersion. We 
find that for t' /t — —0.1, the fiat part of the dispersion 



orthogonal to the antinode remains in close proximity to 
the Fermi level for a larger range of fillings (0.83 < n < 
0.87) when compared with the t' = result. We find 
that the resistivity displays a linear-T dependence, and 
the self energy displays MFL characteristics in the same 
range of doping. We also note that the single-particle 
spectra is less coherent in the center of this doping region 
than it is at the quantum critical filling when t' — 0. 
These observations are consistent with an increase in the 
doping region of marginal Fermi liquid character when 
t'/t < 0. 

There are several different scenarios which may explain 
this behavior. Since the doping region where the van 
Hove singularity is near the Fermi level increases with de- 
creasing t'/t, perhaps the simplest understanding of this 
behavior is that the van Hove singularity pinning gives 
rise to the marginal Fermi liquid behavior—^—. Another 
possibility is that both T and negative t' are relevant 
variables, or variables that, when finite, move the system 
away from a critical point. In addition, they are roughly 
similar in their effect, in that when t' = T = the dop- 
ing region of marginal Fermi liquid character shrinks to 
a point, but the region increases with either increasing T 
or decreasing Thus when t' < 0, the quantum criti- 
cal point may be viewed as moving to negative temper- 
atures so that the quantum critical region broadens to 
allow the linear-T resistivity and the pinning of the fiat 
ky dispersion to the Fermi level to exist for a wider range 
of fillings (compare e.g. the third panel of Fig. [T7] and 
inset of Fig. [T|). However, this interpretation is not con- 
sistent with scaling theory where we expect a finite low 
temperature scale like Tp or T* to emerge for any set of 
parameters that take us away from the quantum critical 
point. Another possibility is that the change in Fermi 
surface topology associated with the van Hove singular- 
ity crossing the Fermi level can introduce Lifshitz singu- 
larities in the free energy^"—. This scenario will mean 
a line of zero temperature critical points in the t' — /i 
plane beyond the quantum critical point as the control 
parameter f is decreased from zero to negative values in 
the t' — fj, plane. For positive t' the Lifshitz transition 
may yield a line of first order critical points in t' — /i 
plane which terminates at the quantum critical point for 
t' = 0. To understand our results, the line of Lifshitz 
transitions for negative values of t'/t must yield an as- 
sociated V-shaped region of quantum criticality which 
becomes flatter as t' decreases as shown on the right in 
Fig. 1171 Finally, another mechanism enhanced by the 
van Hove singularity is the Pomeranchuk instability^^ 
of the Fermi surface where the Fermi surface is distorted 
to break the C4 symmetry of the square lattice. The 
possibility of Lifshitz and Pomeranchuk transitions are 
being studied currently and will be published elsewhere. 

There is some experimental evidence^i in the cuprates 
that there is a change in Fermi surface topology and an 
associated Fermi level crossing of the van Hove singular- 
ity at a doping that is larger than the doping at which 
Tc is maximum, while still being within the dome. On 
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the other hand, we find that the van Hove singularity 
crosses the Fermi level at a slightly smaller doping than 
the optimal doping. This disagreement can be due to 
the other effects present in the real systems e.g. phonons 
that are not incorporated in this model calculation, or the 
strong doping dependence of the strength of the pairing 
interaction^ seen in the simulations. 

The transport provides some additional evidence for 
the van Hove singularity. In our calculations, the low 
energy particle-hole symmetry and the change in sign 
of the thermopower with doping near the critical value 
are both due to the crossing of the van Hove singular- 
ity. However, the doping associated with the van Hove 
crossing differs from that where the thermopower is zero 
due to the anisotropy of the group velocity on the Fermi 
surface. 



V. CONCLUSION 

We explore the role of the van Hove singularity in the 
quantum criticality observed at finite doping in the Hub- 
bard modelj^ Near the quantum critical filling, wc find 
a van Hove singularity due to a flattening of the dis- 
persion near the Fermi level. The motion of the fiat 
part of the dispersion along the antinodal direction is 
anisotropic. The part along the antinode moves contin- 
uously through the Fermi level. The part orthogonal to 
this direction is pinned near the Fermi level at a filling 
where the self energy, transport, and energies^ also dis- 
play marginal Fermi liquid behavior, and the quasipar- 
ticle fraction vanishes i^. Many authors have proposed 
that the van Hove singularity near the Fermi level will en- 
hance superconductivity by enhancing the divergence of 
the bare pairing bubble. Indeed we found previously that 



the superconducting dome surrounds the critical doping 
where the real part of the pairing bubble diverges al- 
gebraically, replacing the Fermi liquid log divergence.— 
However, a simple non-interacting picture with the van 
Hove singularity at the Fermi level doesn't explain the 
quantum critical scaling of the bare dynamic pairing sus- 
ceptibility. We also found previously that a positive t' is 
the control parameter for a first order phase separation 
transition which is terminated by a second order critical 
point. As t' — > this second order terminus is driven to 
zero temperature yielding the quantum critical point.— 
Here we explore the effect of a negative t' , and find that it 
is a relevant variable which increases the extent in doping 
(and chemical potential) of the quantum critical region 
characterized by marginal Fermi liquid behavior. 
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